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SUMMARY 

An experiment has been developed to simulate the wake dynamics of a tur- 
bine rotor stator combination. The experiment consisted of a spoked wheel 
rotating In annular flow, generating simulated rotor wakes. The wakes Impinged 
on a simulated stator, composed of a circular cylinder with a splitter plate 
behind It. The test rig properly simulated the wake shapes and velocity 
dynamics of a rotating turbine stage. The focus of the research was spanwlse 
averaged circumferentially local heat transfer In the circular cylindrical 
leading edge region of a turbine airfoil. The range of Reynolds number based 
on the stator cylinder diameter was 35 000 to 175 000. The reduced blade pas 
sing frequencies, or Strouhal numbers, also based on the stator cylinder dlam 
eter, ranged from 0.63 to 2.50. The wakes were generated by two different 
sets of circular cylindrical bars, 1.59 and 3.18 mm In diameter. Ihe rotor 
could be rotated either clockwise or counterclockwise, causing the wakes to 
Impinge on either the leeward or the windward side of the heated cylinder. In 
addition, grid turbulence was Introduced upstream yielding freestream turbu 
lence of 1.0 to 2.5 percent at the stator location. 

The data presented in this paper, therefore, represent an extensive body 
of local heat transfer coefficients, which can be used to develop models In 
the very Important leading edge region of a turbine airfoil. Ihe wakes pro 
duced only a modest Increase In heat transfer over the baseline no rotor data. 
Ihe data In the presence of the rotor wakes showed an asymmetry from the lee- 
ward to the windward side. Windward side levels of from 30 to 40 percent 
higher than corresponding leeward-side levels were not uncommon, particularly 
at the combination of high Strouhal number and larger wake generating pins. 

Ihe stagnation gauge data were examined against the luVRe parameter over a 
wide range from 8 to 40, and found not to be fully correlated by this often 
used parameter. There were distinct trends with both free stream turbulence 
level and wake generating bar diameter. The results strongly indicate a need 
to separate the wake influence from turbulence effects. 


1NTR0DUC1 LON 

Considerable effort is being directed these days at heat transfer in 
unsteady flow fields. Probably among the oldest unsteady effects studied are 
those which occur naturally In the flow, such as one body In the wake of 


*A shortened version of this report will be presented at the 8th 
International Heat Transfer Conference, San Francisco, California, 
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another or the effect of natural Instabilities. Another area Involves oscil- 
lations of the entire flow field, such as might occur downstream of a pump. A 
third type of unsteadiness, the type of Interest herein. Is a forced, periodic 
and localized, disturbance of the flow, such as occurs where a moving object 
cuts across the flow path, as In a turbine. 

Research on heat transfer In these flow fields Is just beginning. Work 
Includes heat transfer In turbines (refs. 1 to 4), as well as in the wakes of 
spoked wheel or squirrel cage disturbance generators (refs. 5 to 9). The exact 
nature of the moving wakes on heat transfer awaits more research; however, in 
general It can be said that these disturbances enhance heat transfer. The 
passing of the wakes can produce some dramatic local unsteady heat transfer 
effects, as shown by Doorly and Oldfield (ref. 6), Ashworth et al. t (ref. 7). 
The effect on steady heat transfer, especially when combined with normal turbu 
lence, Is not uniformly clear, as shown by Dunn (refs. 1 to 3), for example. 

There Is a growing opinion which suggests that It will be necessary to 
sort out the periodic unsteadiness from the random turbulence. Jayaramon, 
Parlkh, and Reynolds (ref. 10) propose a triple decomposition of the velocity 
vectors. Adamczyk (ref. 11) also recommends changes In turbomachinery closure 
modeling. The extend and detail of what must be done will be evolving for 
sometime. 

The experiment reported herein simulated a turbine rotor wake with a 
rotating spoked wheel. The leading edge of the airfoil In the wake was simu- 
lated with a circular cylinder with a splitter plate behind it. Spanwise 
averaged, circumferentially local heat transfer was measured for angles from 
0° to 79° from stagnation on both the windand leeward side relative to the 
convected wake. Reference data were acquired at various turbulence levels 
with the rotor removed. An attempt Is made to sort out the two effects. 


NOMtNCLATURt 

A cross sectional area, m 2 

A s surface area, m 2 

d diamter, m 

f frequency, Hz 

1 current, amps 

k thermal conductivity, w/m 2 K 
m mass flux, kg/m 2 sec 
Nu Nusselt number, dimensionless 

Pr Prandtl number, dimensionless 

Q heat, power, w 

Re Reynolds number, dimensionless 
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S Strouhal number, dimensionless 
T temperature, K 

Tu unsteadiness parameter, dimensionless 

u average velocity, m/sec 

u' fluctuating velocity, m/sec 
V voltage, volts 

y viscosity, N-sec/m^ 

Subscripts : 
avg average 

aw adiabatic wall 

ref reference 

o stagnation condition 

s static condition 

w wall condition 


R010R WAKE SIMULATION 

lhe simulation employed In this experiment Is Illustrated In figure 1. 

In part (a) a typical turbine stage with the wake- related velocity triangles 
Is shown. Only two points are Important. First, whether one Is looking at 
the rotor passing Into a stator wake or a rotor wake being carried Into a sta- 
tor passage, the effect Is the same. The suction surface sees the wake first 
and the change in Incidence is towards negative Incidence. The simulation 
produced by a circular rod passing across the flow field is shown in part (b). 
The major features are clearly the same. 

The simulation used in the present experiment is illustrated in figure 2. 
A spoked wheel rotates In an annular flow passage producing the desired wakes. 
Since the focus Is on the leading edge region, the downstream airfoil is evpn 
further simplified by using a circular cylinder with a splitter plate behind 
it to eliminate alternate shedding vortices. 

Typical Instantaneous and ensemble- averaged wakes produced In the simula- 
tion rig (more Information on them appears in ref. 5) are shown in figure 3. 
The combination of the turbulence and the periodic unsteadiness Is apparent 
In figure 3(a). The random components disappear in the averaging process, 
figure 3(b). These wakes were cpmpared to compressor wakes measured by Evans 
(ref. 12) and found to exhibit all the major features of depth, width and so 
on. 
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DESCRIPTION OF EXPERIMENT 


The key element In the present experiment Is the segmented copper heat 
transfer cylinder, shown In figure 4. The test section was a 2.62 cm diameter 
circular cylinder composed of 15 Individually controlled heater elements, which 
served as heat transfer gauges. Each heater was made of a solid copper bar 
with a thermocouple (Type E) burled near the surface and a heater element at 
the back. The temperature signal from each element was fed Into an electronic 
control circuit, described In reference 13, which controlled the power to each 
heater such that the temperature of all elements was held constant to within 
0.1 K. Thus, the experiment Is an Isothermal wall boundary condition. Twelve 
of the heater elements were spaced around the cylinder circumference at noml 
nally 11° Increments, as can be seen In figure 4. The two end elements served 
as guard heaters, leaving 10 active heat transfer elements. Both of the span 
wise ends had a guard heater, as shown. The back of the whole assembly was 
also guarded. This element, of course, cannot be seen In the photograph. 

Each of the active heat transfer elements was 2.50 mm wide by 33.0 mm long. 

Each element had a mylar and epoxy Insulating space between them. It can be 
seen from the photograph that the Installation was not perfectly equlspaced. 

(It must be remembered that all 15 elements with heater lead wires and thermo 
couples for each were stuffed Into a package only 2.62 cm In. dlam.) The 
actual center- to-center spacing was measured and these actual angles are used 
In the data reporting. The photograph, figure 4, was taken after the experi- 
ment and shows a small nick In one of the heaters. This occurred midway 
through the experiment, when the smaller pins broke. It appeared to have no 
effect on the subsequent data. 

Ihe circumferentially segmented heated test cylinder was Installed In the 
rotor-wake rig, shown In figure 5, at a distance 2.5 cm downstream of the 
spoked rotor. Thus, the wakes were being generated one stator diameter 
upstream of the stator stagnation line. The flow facility and Its characteri 
sties are described In some detail In reference 5. Only those details essen 
tlal to the present experiment will be described herein. Ihe flow facility 
was at room temperature and pressure annular flow wind tunnel. Ihe flow annu 
lus was 40.6 cm o.d. by 27.1 cm 1 . d . , providing a 0.722 m 2 flow cross section. 
Ihe heated cylinder was one of eight stators equlspaced circumferentially 
around the annulus, each with a 1 .6 mm thick splitter plate behind, as Illu- 
strated In figure 2. The heated stators provided a 20 percent flow blockage. 
The wake generating rotor was a spoked wheel of 24 cylindrical pins, as Illu- 
strated In figure 2. There were two pin diameters, 1.59 mm and 3.18 mm. Ihe 
rotor was driven with a motor and pulley arrangement. It could be driven both 
forward and reverse and could be controlled at speeds from 500 to 7300 rpm. 

The turbulence was controlled with square pattern woven wire screens with 
approximately two thirds open area. Two grid sizes, 10-mesh and 3 mesh, were 
used. (The x mesh terminology Is an American wire manufacturer designation 
standing for x wires per inch and does not translate well to SI units. For 
convenience this jargon will be used herein.) The 3-mesh grid had 1.6 mm wires 
spaced about 8.5 mm apart and the 10 mesh grid had 0.25 mm wires spaced about 
2.5 mm apart. Except for one dataset the grids were located 21.0 cm upstream 
of the stator leading edge. (In one case the 3-mesh grid was located 10.8 cm 
upstream. ) 

As Indicated above, the heat transfer from each element was determined by 
measuring the power delivered to each Individual backside heater and the temp 
erature of the copper element with a thermocouple buried in It. Ihe tunnel 



mass flow rate was measured with a sharp edged orlflced designed to ASME stan- 
dards. This measurement was backed up by treating the Inlet contraction as on 
Isentroplc nozzle. Velocity was measured at two locations, one half way between 
the turbulence grid and the stator plane and the other at the stator stagnation 
plane, using a pilot-static probe. Velocities and turbulence Intensities were 
also measured at the same locations with single-wire hot wire anenometers. 

The rotor speed was measured with a fiber optic probe sighted on the tip of 
one of the rotor pins painted white. This was backed up by a 60- tooth gear 
and magnetic pickup on the rotor drive shaft. 

The data are presented In terms of Nusselt, Reynolds and Strouhal numbers 
and also an unsteadiness or turbulence parameter, all defined below. 

Ihe Reynolds number Is based on a mass- averaged velocity and free-stream 
properties. 
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md 
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The Strouhal or reduced frequency Is based on the rotor circular frequency and 
the same mass- averaged velocity. Expressed In terms of the rotor pin (blade) 
passing frequency It Is: 
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The Nusselt number Is defined In the following manner. 
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The total heat flux Is based on the total power dissipated, VI. The heat loss 
correction has two components; the radiation loss, which was usually quite small, 
and the loss through the Insulating strips, which was occasionally large and 
varied from strip to strip. The Insulation loss was determined analytically, 
based on a model presented In reference 13. The adiabatic wall temperature 
assumed a laminar flow Prandtl number relationship 

T = T - VPr (1 - 1 ) 

aw o’ o s 

and the thermal conductivity was calculated on the basis of the Eckert (ref. 14) 
reference temperature 

T ref ' °- 50 T «) avg * 028 ’s " °' 22 'aw 
Finally, an unsteadiness parameter Is defined as 


u 
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The unsteadiness Is simply the root-mean-square of the velocity flucatlons 
divided by the Integrated time average of the velocity, as measured by the hot 
wires. In the absence of the rotor wakes It is the familiar turbulence Inten- 
sity, Tu. In the wake of the rotor it Includes both the random turbulence and 

the periodic unsteadiness of the wakes, as shown In figure 3. 

An uncertainty analysis, based on the method of Kline and McClIntok 
(ref. 15), was Included In the data reduction program and every dataset 
Included uncertainty estimates. These will be shown on the results figures, 
as appropriate. While the uncertainty varied from run to run, and in the case 

of Nusselt number from heater strip to heater strip, the average values were: 

Nusselt number 6 percent, Reynolds number 1 percent, Strouhal number 1 percent, 
and unsteadiness parameter 1 percent. 

Finally, before initiating the experiment, flow angle measurements were 
made behind the rotor. An example for the 1.59 mm pins is available in 
reference 5. For both pin sizes the combinations of Reynolds and Strouhal 
selected were such that the flow angle was always less than 3°. 


RESULTS 

No Rotor Reference Cases 

In order to establish a reference database, and to compare existing data 
to the literature a series of cases were run with the spoked rotor removed. 

In the first case the heated cylinder was rotated to various angular positions, 
causing different elements to be the stagnation gauge, in order to see if there 
was any gauge bias. The results, shown in figure 6, show very little sensi- 
tivity to gauge position except for one case at the lowest Reynolds number. 

The variation from gauge to gauge was well within the average uncertainty band 
of the stagnation gauge, which is included on the figure. One will note that, 
except at the stagnation point, the data do not line up with the angles. This 
is because, as explained earlier, the gauges are not spaced at equal angles; 
and, when the cylinder Is rotated, the angle relative to stagnation changes. 
There appears to be a slight dip in heat transfer in the stagnation region. 
Since the dip Is within the uncertainty band, one must be cautious in drawing 
conclusions. However, the dip appears to be Independent of gauge selection 
and has been observed by others (ref. 16). 

In order to obtain data over a large angle range (0° to 79°) gauge 8 was 
selected to be the stagnation gauge. It was selected over number three, 
because the uncertainty was always lower. With this choice all the remaining 
data were acquired with gauge 8 on the stagnation line. This left two active 
gauges to the opposite side as a check. In the following figures only the 
eight gauges to one side of stagnation will be shown to avoid confusion. 

The no rotor data over the full Reynolds number range of the experiment 
are shown in figure 7. The error band for each gauge is shown for one of the 
Reynolds numbers, as typical. The data represent two turbulence levels, 
approximately 1 percent and 2 percent. For each level, higher Reynolds numbers 
yielded slightly lower turbulence Intensity. The doubling of turbulence inten- 
sity from 1 to 2 percent leads to a 5 to 10 percent Increase In heat transfer, 
this Increase being greater at the higher Reynolds number. The Increase is 
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pretty uniform for all angles. These data form a baseline reference for the 
rotor-wake cases. 

Despite the large body of data on heat transfer to cylinders In crossflow, 
when one wishes to make comparisons, It Is difficult to find datasets which 
matchup all conditions. For comparison to the present data the experiments of 
Kestln and Wood (ref. 17), VanFossen (ref. 18) and Gledt (ref. 18) were used. 
The comparisons are shown In figure 8, for Reynolds numbers In the range of 
95 000 to 100 000. A third dataset, with higher turbulence level, was added 
to the data already shown on figure 7 for Re = 97 200. These data were 
obtained with the 3-mesh grid only 10.8 cm from the heater cylinder, as was 
mentioned above. Although the data of Kestln and Wood (ref. 17) were actually 
mass transfer sublimation data, the flow conditions most closely match the 
present data. The comparison Is quite good except at the large angles. The 
data of Van Fossen (ref. 18) agree very well but his data only go out to 50°. 
The turbulence for these data were generated with parallel wires rather than a 
square mesh grid. The comparison with Gledt (ref. 19) Is less satisfactory 
both In level and angular dropoff. Gledt's data were uniform heat flux rather 
than uniform temperature. Papell (ref. 20) showed that uniform heat flux data 
dropoff more slowly than uniform temperature, which could explain the observed 
angle variation. The level, especially at stagnation Is harder to explain, 
except that Gledt estimated his turbulence Intensity, while all the others 
measured It. In general, within experimental uncertainty, the present data 
agree well with the literature, well enough at least to draw valid comparisons 
between the no-rotor and rotor-wake data In the same experiment. 


Rotor Wake Data 

A very extensive dataset was acquired In the present experiment, more 
than can be presented herein. A selected subset Is shown In figure 9(a) to 
(d). Two midrange Reynolds numbers and all three Strouhal numbers are pre- 
sented. Each plot Is for a different combination of pin diameter and turbu- 
lence grid, four In all. In addition, each plot Includes data from wakes 
produced by both forward and reversed motion of the spoked rotor. A small 
diagram of the geometry of the flow field Is shown on each plot. For the data 
to the left of center the wake Is convected In at the heat transfer gauges 
( 1 . e . , the windward side) and to the right of center the gauges are In the 
shadow of the convected wake ( 1 . e . , the leeward side). Recall, the wake flow 
angles are always less than 3°. In terms of the simulation the left hand 
(windward) data would represent the suction surface of a turbine airfoil and 
the right hand (leeward) the pressure surface. 

Examining figures 9(a) to (d), one Is struck by the strong similarity 
from case to case. The effects of Strouhal number are small but systematic. 
Although the effects are on the same order of the uncertainty bands, the fact 
that they are so consistently systematic Is probably relavant. There Is a 
definite asymmetry In the data, especially at the higher Strouhal numbers and 
with the larger pins. 

In order to get a better appreciation for the various effects, a cross 
plot of the data was used. The data for one Reynolds number (125 400) and one 
Strouhal number (1.27) are plotted for all four pins and grid combinations. 

In addition the no-rotor data are plotted, Including the highest turbulence 
(grid at 10.8 cm) case. To aid In comparison the no- rotor data are plotted to 
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both the right and left of zero on figure 10. They are not two datasets. 

Thus, gauge 7 data appears at ±10°, gauge 6 data at ±25°, and so on. The 
average uncertainty band for the stagnation gauge is shown on the figure. 

Many of the observations to be made fall within this band. 

First, let us focus on the data associated with the 10-mesh grid. In 
the stagnation region the rotor-wake data show a very modest (6 to 7 percent) 
Increase In heat transfer, which Is within the uncertainty. In fact, the 
small-pin data are higher than the large-pin, which does not seem right. 
Nevertheless, the data are higher. On the leeward side these modest increases 
persist over the whole angle range. On the windward side, as the flow pro- 
gresses to higher angles, a more substantial wake-induced Increase In heat 
transfer occurs, reaching almost 50 percent over the no- rotor data at 79°. 

This Is, of course, near separation. However, except at the 79° point, all 
of the 10-mesh data fall below the no-rotor 3-mesh data sets. This, despite 
the fact that the wake related turbulence/unsteadiness rms Is approximately 
5 percent for the 1.59 mm pins and 9 percent for the 3.18 mm pins. In fact, 
an almost doubling of the rms unsteadiness seems to have almost no effect on 
average heat transfer. This Important result points out the danger of treating 
the wake-induced unsteadiness as a simple turbulence intensity. 

The Insertion of the 3-mesh grid at the 21.0 cm upstream location produced 
about a 70 percent Increase In the baseline turbulence. This Increase shifted 
all of the spanwlse-averaged heat transfer data, both no-rotor and rotor wake, 
up uniformly, about 13 percent at stagnation. Curiously, Increasing the base- 
line turbulence appears to have little or no effect on the rms turbulence/ 
unsteadiness of the rotor wakes, but It does affect the heat transfer. Finally, 
one notes that Increasing the grid turbulence to a modest level of 2.84 percent 
produces a higher heat transfer level than all the rest of the heat transfer 
data. 

It has become common to express cylinder In crossflow heat transfer data 
as Nu/-s/Re verses Tu\/Re". Kestln and Wood (ref. 17) and Lowery and Vachon 
(ref. 21) proposed quadratic correlations for the stagnation line heat transfer 
which are often used. Lowery and Vachon (ref. 21) also showed that the same 
trends as at stagnation existed for overall heat transfer, but, of course, at 
lower average values. In reference 5 the present authors suggested that over- 
all heat transfer to a circular cylinder In the wake of a spoked rotor could 
be correlated In the same manner by treating the rms flucuatlons behind the 
rotor as turbulence. Following up on that suggestion, the stagnation gauge 
data for all cases are shown in these variables on figure 11 along with the 
data of Lowery and Vachon. The no rotor data agree well with data of Lowery 
and Vachon. The rotor-wake data do not agree and they do not correlate in 
these variables. This Is not too surprising In light of the previous figure. 
What was thought to be a correlating trend In reference 5, was not. The dif- 
ference In the present experimental conditions from the previous Is the addl 
tlon of 3.18 mm pins. They were not used In reference 5. Since the larger 
pins double the rms flucatlons but make very little effect on the heat trans- 
fer, Tu' cannot be a correlating parameter. Actually, if one were to Ignore 
the wake and use the grid turbulence level, the data would probably fall within 
the wide band of the no rotor data, shown In figure 11. This representation 
would be wrong, however, since the wakes do have an effect on stagnation heat 
transfer, albeit small. 
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The main effect of the wakes on the spanwise average heat transfer Is the 
angular asymmetry. This effect Is best shown In figure 12, where the ratio of 
wlndward-to-leeward heat transfer Is plotted against angle from stagnation for 
the data of figure 9. There Is quite a bit of scatter In this data, neverthe- 
less clear trends do exist. Lines were drawn roughly through the data to aid 
In Interpretation. In this representation It can be seen clearly that pin 
size and Strouhal number (rotor speed) make a difference. There appeared to 
be no special Influence of grid turbulence or Reynolds number on this asym- 
metry. Although there Is a definite trend at all angles, at angles below about 
40° the asymmetry was less than 5 percent, which Is within the uncertainty.- 
At the lowest Strouhal number, 0.62, the asymmetry ranged from 5 percent with 
the small pins to 10 percent with the large pins. At the highest Strouhal 
number, 2.51, the asymmetry range was 35 to 50 percent at the higher angles. 

It should be noted that on an airfoil the circular leading edge Is faired 
smoothly to the rest of the airfoil and thus separation effects are minimized. 
The large effects seen herein at large angles should not be present on an air- 
foil but the asymmetry will be present. Further Investigation on airfoils Is 
warranted . 


CONCLUSIONS 


An experiment was performed In which the wake effects of a turbine stage 
were simulated using a spoked-wheel rotor and downstream circular cylinder as 
a stator leading edge. Spanwise average circumferentially local heat transfer 
coefficients were measured over a range of Reynolds and Strouhal numbers and 
at two levels of grid-generated turbulence. The rotor could be rotated In 
either direction, resulting In the wakes being convected both Into and away 
from the heated surface. The primary conclusions from this experiment are as 
follows: 

1. The wakes Induced by the spoked rotor produce a modest Increase In 
(heat transfer over the no-rotor reference case. 


2. Increasing the baseline turbulence shifts both the no rotor reference 
case and the rotor-wake heat transfer data up uniformly. 


3. The wake- generated rms unsteadiness Is not a correlating parameter 
and should not be classified as turbulence. The data strongly support the 
need to separate out the effects of random turbulence and periodic unsteadiness 
and treat them separately. 


4. The convected wakes produce a significant asymmetry In heat transfer 
around the leading edge of the heated cylinder from the side Impacted by the 
wakes to the side shielded from the wakes. This asymmetry Increases with 
Strouhal number and wake size. 
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(b) Simulation. 

Figure 1. - Simulation of turbine wakes using a spoked wheel disturbance 
generator. 



Figure 2. - Concept of rotor wake heat transfer experiment. 
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(b) Ensemble averaged wake. 

Figure 3.- Typical wakes generated by a spoked wheel disturbance 
generator. 
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Figure 4. - Photograph of circumferentially segmented copper heater. 



Figure 5. - Rig flow schematic. 
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Figure 6. - Effect of circumferential position of stagnation 
gage on Nusselt number distribution. 
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Figure 7. - Spanwise average circumferential Nusselt number 
distribution as a function of Reynolds number and turbu- 
lence level with the rotor removed. 
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Figure 8. - Comparison of present data (rotor removed) 
with other experiments. 
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(b) 1. 59 mm Pins and 3 mesh grid. 

Figure 9. - Effect rotor pin wakes on spanwise average 
circumferential heat transfer. 
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Figure 10. - Effects of various types of unsteadiness on the 
circumferential distribution of heat transfer. 
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Figure 11. - Comparison of rotor wake stagnation heat transfer data to frequently used correlation. 




Figure 12. - Heat transfer distribution asymmetry resulting 
from wake size and rotational speed. 
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